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Abstract 

An advanced engineering computational model has 
been developed to aid in the analysis of chemical rocket 
engines. The complete multispecies, chemically reacting 
and diffusing Navier-Stokes equations are modelled, 
including the Soret thermal diffusion and Dufour energy 
transfer terms. Demonstration cases are presented for a 
1030:1 area ratio nozzle, a 25 lbf film -cooled nozzle, and a 
transpiration -cooled plug-and-spool rocket engine. The 
results indicate that the thrust coefficient predictions of the 
1030:1 nozzle and the film-cooled nozzle are within 0.2 to 
0.5 percent, respectively, of experimental measurements. 
Further, the model’s predictions agree very well with the 
heat transfer measurements made in all of the nozzle test 
cases. It is demonstrated that thermal diffusion has a sig- 
nificant effect on the predicted mass fraction of hydrogen 
along the wall of the nozzle and was shown to represent a 
si g nificant fraction of the diffusion fluxes occurring in the 
transpiration-cooled rocket engine. 

Nomenclature 

Ordinary mass diffusion coefficient 
Nozzle vacuum thrust coefficient 
Characteristic velocity, ft/sec 

Thermal diffusion coefficient of species i 
Enthalpy of species i 
Diffusion mass flux of species i 
Turbulent mass flux of species i 
Nozzle specific impulse, lbf-sec/lbm 
Effective thermal conductivity 
Heat flux 

Fluid static pressure 
Nozzle chamber stagnation pressure 
Fluid temperature 
Turbulent diffusion coefficient 

Dufour energy transfer coefficient 
Mass fraction of species i 

* Member AIAA 

♦♦Professor, Mechanical Engineering Department 
Associate Fellow, AIAA 


Introduction 

The computation of rocket-nozzle performance is of 
critical importance to the rocket-engine industry. Often, 
the rocket engine manufacturer needs to compare various 
engine operating systems before deciding on a specific 
configuration. Two items that are of critical importance are 
nozzle performance and nozzle heat transfer characteris- 
tics. Often these two parameters are interrelated. This is 
most apparent in cases involving transpiration cooling or 
film cooling. Under these circumstances, the actively 
cooled boundary layer serves to thermally (and perhaps 
chemically) protect the nozzle wall from the exhaust gas- 
ses of the core flowfield. In such conditions, the coolant 
allows the rocket engine to operate at conditions otherwise 
not attainable, such as higher chamber pressures or lower 
thrust levels. The engine designer must then ascertain 
whether the increased operating pressure or lower thrust 
level is worth the decreased performance of the engine 
because of active cooling. Consequently, a need exists for 
a comprehensive rocket nozzle prediction capability that 
can accurately determine performance and heat transfer 
characteristics for design and trade-off studies. 

A different need exists for individuals interested in 
improving rocket nozzle performance predictions and 
modelling. Here, it is important for the individual to ascer- 
tain the primary limitations of the prediction technique. A 
model which would include as many of the potential per- 
formance and heat transfer effects as possible and also 
allow the examination of individual effects, such as spe- 
cies diffusion or chemical kinetics, would be most desir- 
able, with the additional constraint that the solution should 
be reasonably insensitive to the computational grid 
employed for the analysis. 

Present Study 

Modelling of the multispecies rocket-nozzle flowfield 
equations was accomplished by modifying the compress- 
ible Navier-Stokes code, Proteus. 1 ' 3 The Proteus code 
included the ability to model the time-averaged Navier- 
Stokes equations for a single-species perfect gas (constant 
Y). For a fluid with the properties of air, the Proteus code 
also included the ability to model a temperature-dependant 
specific heat, thermal conductivity, and viscosity. 

In the present study, inclusion of the multiple species 
advection, production (i.e., chemical kinetics), and diffu- 
sion terms, with full coupling in the total energy equation, 
was developed. Further, species and heat transfer pro- 
cesses by diff usion were modelled with consideration 
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given to the interrelation between species and heat transfer 
mechanisms (i.e., the Soret and Dufour effects). 

The code was further enhanced by the inclusion of the 
following boundary conditions: 

• Characteristic equations at the boundaries; 

• Mixed subsonic-supersonic exit conditions; 

• Coupled temperature/mass flow wall conditions. 

The resulting fluid-flow model that incorporates all of the 
above modifications to the Proteus code was called Tethys. 

Governing Equations 

The multispecies Navier-Stokes equations governing 
the flow of a fluid are similar to their single-species coun- 
terpart, with the exception that additional species-transfer 
equations must be solved. The species-transfer equations 
contain advective, source (due to chemical reactions), and 
diffusive components. The governing equations can be 
found in Reference 4 and, for brevity, are not repeated 
here. A unique aspect of multispecies fluid flow is that the 
constitutive equations governing the flow of heat and spe- 
cies transfer are significantly more complicated than those 
of a single-species gas. 

The constitutive equations of heat transfer are compli- 
cated by the presence of multi-species interactions as heat 
can transfer diffusively by means of thermal conduction, 
interdiffusion, and Dufour heat transfer, as shown below: 

q = -kVT + £h i o i +jl)+i 

i = l i=lj = l L i J - 1 

Conduction Interdiffusion Dufour heat transfer 

The transport of species can be the result of ordinary 
diffusion and thermal (i.e., Soret) diffusion, with the addi- 
tional complication that the Boussinesq approximation for 
turbulent mass transfer in not completely analogous to its 
laminar counterpart, as illustrated in the following equa- 
tion: 


n-l 

j = l 

Ordinary Thermal Turbulent 

diffusion diffusion diffusion 


A transformed cylindrical coordinate system was used 
for accurate boundary description and to vary the resolu- 
tion of the computational grid. The specific details of the 
transformation technique employed can be found in Refer- 
ence 5, with the unique details of the transformation 
method outlined as follows: 

• All radius terms contained within axial derivatives are 
evaluated as spacial averages; 

• The differencing of the metrics are never beyond the 


boundary at which a dependent variable has been 
evaluated; 

* The pressure gradient in the radial momentum equa- 
tion is described in weak conservation form. 

The reformulation using the above guidelines results 
in discretized expressions that are both more robust and 
more accurate, based on fundamental flow situations, than 
the method used in the Proteus code. 


The application of the reacting, multispecies diffus- 
ing, discretized flowfield equations will be shown for three 
unique hydrogen/oxygen rocket engines that were tested at 
the NASA LeRC rocket-engine test facilities. First, a high 
area ratio (1030:1) rocket nozzle was analyzed and dis- 
cussed in detail. The application of such a nozzle has been 
envisioned to include orbital transfer, lunar, and planetary 
missions. An extensive experimental database exists that 
includes performance, heat flux, and static wall pressure 
measurements. A second hydrogen-oxygen rocket case 
that will be presented is a film-cooled, low thrust (25 lbf) 
engine that has been proposed for satellite thruster applica- 
tions and for the orbital maintenance of a man-tended 
space station. The uniqueness of this engine concept was 
the extensive film cooling with hydrogen and operation at 
an overall stoichiometric mixture ratio while having an 
oxygen/hydrogen mixture ratio of 20 in the core. Static 
wall pressures and global heat transfer and engine perfor- 
mance data are available for comparative purposes. The 
final configuration examined is a transpiration-cooled 
plug-and-spool rocket engine. The uniqueness erf this con- 
figuration was that a discrete section of the rocket has tran- 
spiration hydrogen coolant along the wall. 


A high area rath 
ure 1) was tested at 
Gaseous hydrogen and oxygen were injected into a com- 
bustion chamber that had a subsonic area ratio of 4.22:1. 
In the chamber, the propellants were burned to near com- 
pletion, with the combustion efficiencies in each test rang- 
ing from 93 to 99 percent 


> hydrogen/oxygen rocket nozzle (Fig- 
the NASA Lewis Research Center. 8 



Fig. 1. LeRC 1030:1 nozzle dimensions. 


Four different nozzle cases, with various assumptions. 
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were considered and are shown in the following table. The 
analysis cases varied from a chemically frozen, inviscid 
case, to a chemically reacting, multi-species diffusion 
case. 


Table 1. Analysis methods used for the 1030: 1 nozzle. 


Analysis 

Case 1 

Case 2 

Case 3 

Case 4 

Chemical kinetics 

No 

No 

Yes 

Yes 

Viscous flow 

No 

Yes 

Yes 

Yes 

Ordinary diffusion 

No 

No 

Yes 

Yes 

Thermal diffusion 

No 

No 

No 

Yes 

Interdiffusion 
heat transfer 

No 

No 

Yes 

Yes 

Dufour transfer 

No 

No 

No 

Yes 


For analytical comparison, several important assump- 
tions were made concerning the environment under which 
the rocket nozzle performed. The major assumptions 
regarding the flowfield included complete mixing at the 
nozzle inlet, the magnitude of the ambient pressure at the 
nozzle exit, and laminar flow throughout the nozzle. 

The flow at the inlet of the nozzle section was 
assumed to be fully mixed and in chemical equilibrium, 
leaving the lack of complete mixing as a ‘correction fac- 
tor’ to be applied to the engine performance measure- 
ments. 

The ambient pressure at the nozzle exit can be an 
important parameter as the nozzle could experience some 
level of ‘boundary layer feedback’ where the subsonic 
flow region of the boundary layer is affected by the pres- 
sure at which the nozzle exhausts to. In turn, this can affect 
the pressure on the nozzle wall for some distance upstream 
of the nozzle exit plane. Since the nearest static pressure 
tap on the nozzle was located 2.7 inches from the exit, a 
q uail extrapolation (4.1 %), using the interior nozzle static 
pressures was required to establish the nozzle ambient 
pressure. It should be noted that the ambient pressure 
applied in this analysis was 30 % lower than the capsule 
pressure, thus indicating that the nozzle was somewhat 
overexpanded. 

The nozzle flow was assumed to be laminar through- 
out the flowfield. Extensive studies and correlation with 
experimental and analytical predictions indicate that the 
rapid acceleration of the exhaust gasses in the throat 
region of the nozzle ‘laminarizes’ the flowfield. 7, 9 Addi- 
tional studies have indicated that the laminarized profile 
was stable throughout the supersonic region of the noz- 
zle. 10 The transonic and supersonic region of the nozzle 
was the area of greatest concern, as this was the region 
where the growth of the boundary layer was rapid and the 
divergent region of the nozzle was also the location where 


all of the experimental wall measurements were made. 
Consequently, the assumption of laminar flow throughout 
the nozzle was reasonably justified. 

Boundary conditions employed in the analysis used 
parameters measured in the experimental program such as 
the chamber pressure and chamber mixture ratio, and the 
wall temperature. The nozzle ambient pressure was used 
as a boundary condition when the normalized contravari- 
ent velocity was less than the speed of sound. At the noz- 
zle inlet, the flow was assumed to be in chemical 
equilibrium. 

Constant reference-plane characteristic relations were 
used at all boundaries. At the inlet and exit, the character- 
istic relations included viscous, chemically reacting terms. 
Along the nozzle centerline and wall, a chemically-frozen, 
inviscid characteristic relation was used. When applied to 
a no-slip wall, the characteristic equation reduces (at 
steady state, and in the limit as the grid spacing reduces to 
zero) to the assumption that the normal pressure gradient 
is equal to zero. Additional boundary conditions along the 
centerline (and along the wall in the inviscid analysis) 
included the assumption that the gradient in enthalpy and 
entropy were zero. 

rnmnntational Grid U sed for Analytical Study 

In Figure 2, the 91X51 numerical grid used for the vis- 
cous analysis of the 1030:1 nozzle is shown. A slightly 
more coarse grid was used in the inviscid analysis case but 
was essentially indiscernible to casual visual inspection, 
and is not presented. 



Fig. 2. Computational grid used for LeRC 1030:1 nozzle. 


The most unique aspect of the computational grid 
shown in Figure 2 is that the centerline length is signifi- 
cantly shorter than the nozzle length. This ‘centerline cor- 
ing’ of the computational grid takes advantage of the fact 
that, in supersonic flow, information can only propagate 
within a discrete region (the domain of influence) and 
computer tests run with coarser grids indicated that the 
grid shown in Figure 2 was well downstream of the 
domain of influence of the supersonic core on the nozzle 
wall. The net result of this ‘coring’ in the centerline region 
was that it allowed significantly more grid points to lie in a 
physical location where they would be of use to resolve 
the flowfield and the conditions that exist along the wall. 
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Mach Number Conto ur Predictions 

In Figure 3 to 6, Mach number contours as predicted 
for Case 1 and 4 are shown, with the other two viscous test 
cases having profiles qualitatively similar to those pre- 
sented in Figures 5 and 6. Close inspection of Figure 3 
reveals that the flow near the wall of the nozzle in the 
supersonic region undergoes a large amount of compres- 
sive turning. The result of the compressive turning is that 
the Mach number contours are nearly parallel to the wall. 



Fig. 3. Throat region Mach number predictions. Case 1. 



Fig. 4. Predicted Mach number contours. Case 1. 



Fig. 5. Throat region Mach number predictions. Case 4. 



Fig. 6. Predicted Mach number contours, Case 4. 

It is also observed that the centerline boundaries 
ap pea r to over dampen the Mach number contours in the 


supersonic regions of the flow. This may be because of the 
comp utatio nal grid employed for the analysis as the con- 
stant \ lines go from left -to- right as the physical location 
goes from the center to the wall, as illustrated in Figure 2. 
However, in supersonic flow, the centerline can only be 
influenced by conditions existing within its Mach cone, 
which is in the leftward direction. As a result, any relation 
developed between a centerline point and an interior point 
located outside of the domain of dependance is physically 
unrealistic. In the computational grid employed in this anal- 
ysis, this aspect is far more pronounced than in an analysis 
that employs orthogonal grid structures at a centerline, 
although even an orthogonal grid violates the principle of 
the domain of dependance. Fortunately, the analysis of the 
centerline region appears excellent in the transonic/ slightly 
supersonic region and the error that occurs near the end of 
the centerline region will have little effect on the wall con- 
ditions. 

Wall Static Pressure and Heat F lux Comparison 

In Figure 7 and 8 a comparison for Case 4 of the pre- 
dicted and measured pressure and heat flux distributions in 
the supersonic region of the nozzle (the only region where 
measurements were made) are shown. The predictions and 
measurements are in excellent agreement. 



Fig. 7. Comparison of measured and predicted wall 
static pressures. 

It is especially noteworthy that the ‘kink’ in the slope at 
the last three points of wall pressure measurements shown 
was correctly predicted. 



The heat-flux predictions are slightly higher than 
experimental measurements in the near throat region, fol- 
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lowed by predictions that are slightly lower than experi- 
mental measurements. Excellent agreement of 
experimental and predicted heat fluxes occurs from about 
20 in. axial distance to the exit of the nozzle. 

Prediction of Hydrogen Mass Fraction Along the Wall 

In Figure 9 the mass fraction distribution of diatomic 
hydrogen predicted to occur along the wall of the LeRC 
1030:1 nozzle is presented for all of the viscous nozzle 
analysis cases performed. As expected, the concentration 
of hydrogen along the nozzle wall for the frozen analysis 
(Case 2) was constant throughout the length of the nozzle. 
On the other hand, the mass fraction of hydrogen for the 
chemically-reacting, ordinary-diffusion case (Case 3) 
increases slightly at the inlet due to the recombination of 
monatomic hydrogen, and then decays towards the fully- 
mixed, frozen mass fraction as the fluid undergoes a com- 
bination of species mixing and freezing of the chemical- 
reaction mechanisms. The most pronounced difference in 
the predictions of hydrogen mass fraction occurs in Case 
4, where Soret mass transfer was considered. In this case 
the hydrogen mass fraction immediately drops to one-half 
its inlet value, then slowly increases as the boundary layer 
thickens in the combustion chamber. The hydrogen mass 
fraction reaches a local maximum near the throat where 
the strength of the Soret mass transfer term was weak from 
the presence of a modestly cooled, ceramic coated wall. 
The hydrogen mass fraction then decays to a local mini- 
mum as the Soret mass transfer term increases from the 
rapidly decreasing wall temperatures in this region. 
Finally, as the ceramic-coated section ends and the heat 
fluxes diminish, the hydrogen mass fraction again 
increases as the temperature gradients at the wall decrease 
and the boundary layer thickens, essentially remixing the 
fluid. This partial-remixing process continues as the noz- 
zle temperature gradient continues to decrease throughout 
the extent of the expansion region of the nozzle. 



Fig. 9. Hydrogen wall mass fraction predictions. 
Performance Predictions 

A comparison of predicted performance parameters 
for the four numerical test cases analyzed is shown in 
Table 2. The procedure used to determine nozzle thrust 
and mass flow rate corresponds as nearly possible to that 


recommended in the JANNAF nozzle analysis proce- 
dure. 11 This method involved the numerical integration of 
the pressure forces and axial-momentum flux across the 
throat of the nozzle, which in this case corresponded to the 
first grid line that extended into the divergent portion of 
the nozzle, and summing the pressure and drag forces 
along the divergent surface of the nozzle. 

Table 2. Performance predictions for the 1030:1 nozzle. 



Case 1 

Case 2 

Case 3 

Case 4 


488.2 

481.6 

493.6 

494.1 

C F 

1.959 

1.924 

1.941 

1.942 


Several conclusions can be drawn. The first is that 
chemical kinetics seems to have the largest single effect on 
nozzle performance, with the performance difference 
between the frozen and kinetic cases (Cases 2 and 3, 
respectively) being about 2.5 %, where the frozen case is 
lower. Examination of Cases 1 and 2 reveals that viscous 
effects account for approximately a 1.4 % decrease in noz- 
zle performance. The effects of multidiffusion mass and 
heat transfer, the difference between Case 3 and 4, account 
for a predicted increase in nozzle performance of 0.1 %. 

Comparison of Performance Predictions 

In Table 3, a comparison between the predictions 
made with Tethys, Case 4, and experimental measure- 
ments is shown along with previous performance predic- 
tions using other models, namely TDK and RPLUS. 12,13 
The TDK prediction method involves an inviscid (one- 
dimensional/transonic/MOC)- boundary layer type of 
analysis with chemical kinetics in the inviscid region 
while RPLUS is a chemically reacting Navier-Stokes 
model which was executed using a 270X60 computational 
grid. The measured specific impulse has been corrected for 
incomplete mixing and/or combustion by the use of an 
inviscid, chemical-equilibrium correction factor. 

A comparison of the predictions to the experimental 
measurements indicates that the Tethys flow model is in 
far better agreement with the experimental measurements 
than other computational models. Here it is seen that the 
predicted thrust coefficient and specific impulse from the 
Tethys model are about 2 % closer to experimentally mea- 
sured values than the other nozzle-prediction models. 


Table 3. Perfoimance comparisons for the 1030: 1 nozzle. 



Measured 

Predicted/Measured, % 

Tethys 

TDK 

RPLUS 

C F 

1.946 

-0.21 

-2.7 

-1.4 

hp 

495.6 

-0.30 

-2.8 

-3.7 


The differences between measured and predicted val- 
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ties of thrust coefficient and the differences between the 
measured and predicted specific impulse should be nearly 
the same. The difference between the two values reflects 
the predicted value of the discharge coefficient of the noz- 
zle. In the Tethys model, the predicted discharge coeffi- 
cient is equal to 0.991 while the TDK and RPLUS models 
predict 0.99 and 0.977, respectively. Consequently, the 
proximity to which the discharge coefficient of the nozzle 
is equal to unity lends additional credence to the accuracy 
of the Tethys code. 


25 lbf Film-Co oled Rocket Engine 

A low-pressure (75 psia), low-thrust (25 lbf) hydro- 
gen/oxygen rocket engine was tested at the NASA LeRC 
as part of a nozzle contour evaluation program. 14, 15 A 
schematic of the nozzle geometry is shown in Figures 10 
and 11. In Figure 10, a schematic of the nozzle geometry 
and test apparatus is shown. In Figure 12 a view of the 
film- cooling and injector apparatus is presented. From 
Figure 11, the hydrogen film acts as a wall jet, the hydro- 
gen being injected along the surface of the nozzle, and ini- 
tially separated from the core by a solid- walled sleeve. 



Fig. 10. Film cooled nozzle geometry. 


Sparkplug Fuel splitting washer 

/ Sleeve 



Combustion 

chamber 


Hydrogen film 


Hydrogen coolant slot 
injection 

Oxygen injection 

Fig. 11. Injector and film-cooled sleeve assembly. 


Computational Grid Used for Film-Cooled Nozzle 


The 91X51 computational grid employed for the film- 
cooled nozzle analysis is presented in Figure 12. 

Every other computational line shown 



Axial length, in. 

Fig. 12. 91X51 grid used for film-cooled nozzle study. 


Film-Cooled Nozzle A nalysis Assumptions 

Three critical assumptions in the analysis of the 25 
lbf, film-cooled rocket were made so that analytical results 
could be generated in an expedient manner. The critical 
assumptions were: (1) The flow was turbulent and could 
be reasonably modelled by an algebraic turbulence model, 
(2) the constituency of the core fluid flow of the nozzle 
was known, and (3) the radial and tangential velocity dis- 
tribution at the nozzle inlet was zero for both the film and 
the core flow. 

The flow was assumed to be turbulent in the film/core 
interaction region of the combustion chamber, based on: 
the Reynolds number of the hydrogen coolant in the sleeve 
was about 7000, the Reynolds number of the core flow 
(based on the sleeve thickness) was about 5000, and that 
combustion creates an appreciable level of freestream tur- 
bulence. 

Further down the engine, the flow may have experi- 
enced a significant amount of laminarization, both along 
the wall and in the mixing layer, as the flow was acceler- 
ated through the converging/diverging portion of the noz- 
zle. This flowfield laminarization may not be adequately 
modelled by an algebraic turbulence model. However, 
since the nozzle was short, the overall shear stresses acting 
on the nozzle are relatively minor and contribute only to a 
small degradation of the total nozzle thrust The model 
used in this analysis corresponded to a model applied by 
Gortler to obtain the proper variation of the turbulent vis- 
cosity in a mixing layer, with the near wall analysis being 
performed using the Baldwin-Lomax algebraic turbulence 
model. 16, 17 

For expediency, the assumption of a known chemical, 
thermal, and velocity distribution in the core-flow region 
needed to be made. As this was deemed to be a critical 
assumption in the analysis, two variations of the core flow 
inlet conditions were examined; one which assumed a 
fully-mixed core, and another which was based on temper- 
ature and species measurements of a similar injector/com- 
bustor. 18 

Because the hydrogen film was injected through a 
three-dimensional channel configuration before injection/ 
mixing with the core flow, it is certain that some distribu- 
tion of radial and tangential velocity was present at the 
inlet region of the hydrogen film-cooled section. However, 
these velocity components were assumed to be of minor 
importance compared to the axial component of the veloc- 
ity. Similarly, the oxygen and hydrogen propellants in the 
core region of the combustion chamber were injected radi- 
ally and subjected to the effects of a blunt spark-plug igni- 
tor. Under these conditions, the core flow would 
undoubtedly have some distribution of radial and tangen- 
tial velocity for some distance downstream of the region 
where the propellants are injected Adequate modelling of 
these conditions would require a three-dimensional analy- 
sis which was beyond the capabilities of the Tethys model. 
Instead, the sensitivity of nozzle performance to the 
degree of core mixing was examined. 
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Boundary Conditions for the Film-Cooled Nozzle 

Experimental values of the inlet mass flow rates, wall 
temperature, and exhaust pressure of the nozzle were used 
as boundary conditions in the analysis. To complete the 
boundary conditions, constant reference-plane characteris- 
tic relations were used at the inlet, exit, and nozzle center- 
line, with the exception that in the region around and 
downstream erf the shock wave which strikes the center- 
line. A these boundary points, reflective boundary condi- 
tions were used. 

Two different cases were analyzed, the variation being 
the assumed chemical profile at the inlet region of the core 
flow in the nozzle. One analysis assumed the flow to be 
fully mixed. The other assumed that the flow was strati- 
fied, as presented in Figure 13, and was based on tempera- 
ture and species measurements in a similar injector. 1 



Fig. 13. Oxygen/hydrogen distribution used for 
the stratified core case. 


Mach Number Contour Predictions 

In Figure 14, the predicted Mach number distribution 
for the case of a fully-mixed core flow is presented. Here it 
is observed that a rather weak shock wave coalesces at the 
nozzle centerline and causes the Mach number to decrease 
from a value of around 3 to about 2.5. Also to be observed 
is the distorted nature of the Mach number profile in the 
highly supersonic region of the nozzle, where it is higher 
(at a given axial location) near the surface of the wall than 
the centerline. It is evident that the shock waves had a far 
greater affect near the nozzle centerline region than near 
the wall region. The Mach number distribution in the case 
of stratified flow has nearly the same qualitative properties 
as the fully-mixed case. 



Fig. 14. Mach number contour predictions with the 
fully-mixed core case. 


Predicted Steam Mass Fraction Contours 

In Figure 15 and 16, a contour mapping of the mass 
fraction of steam (H 2 0) is shown for the fully-mixed core 
and the stratified core cases. In examining these two fig- 
ures, numerous flow features become evident One is that 
in both the fully-mixed and the stratified-core cases, fluid 
mixing is essentially non-existent in the supersonic region 
of the nozzle. This is entirely as expected, since the fluid 
convective speeds increase rapidly in the transonic portion 
of the nozzle, while the turbulence quantities tend to ‘lam- 
inarize’. Another important observation is that the strati- 
fied flow-case exhibits a much lower mass fraction of 
steam in the core region of the flow throughout the nozzle. 
Since the creation of steam is associated with a large 
increase in thermal energy, a significantly higher perfor- 
mance can be expected from the fully-mixed core case. As 
will be shown later, this is exactly the situation. 



Fig. 15. Steam mass fraction predictions using the 
fully-mixed core assumption. 



Fig. 16. Steam mass fraction predictions using the 
stratified core assumption. 

Predicted and Experim ental Heat Transfer Rates 

Global measurements of wall heat-transfer rates were 
made. The predicted distribution of wall heat flux is shown 
in Figure 17 for the case of stratified flow, with the fully- 
mixed case having a similar heat-flux profile but possess- 
ing a somewhat larger heat-flux peak. A noteworthy obser- 
vation is that the peak in the nozzle heat flux is 
appreciably upstream of the nozzle throat. This peak heat- 
flux location is well upstream of the location normally 
observed due to nozzle curvature. Based solely on nozzle 
curvature, the peak in heat flux would be expected to occur 
approximately 0.04 in. upstream of the throat. 19 The pro- 
file seems to suggest injector streaking. The streaking 
probably occurs because the momentum of the hot-film 
core was significantly greater than that of the hydrogen- 
film coolant, by about 20 times. The core flow thus has the 
tendency towards direct impingement onto the wall since 
the low momentum of the coolant has little ability to force 
the core flow to turn with the geometry. Consequently, the 
fluid streamlines become concentrated very near the wall 
of the nozzle. 
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Fig. 17. Wall heat fluxes predicted for the film-cooled 
nozzle using the stratified core assumption. 

Another noteworthy aspect of the heat-flux predic- 
tions shown in Figure 17 is the local minimum, followed 
by a local maximum that occurs just downstream of the 
throat. In fact, the sharp increase in heat flux corresponds 
precisely to the region where the geometric transition from 
the circular-are throat to the conical divergence occurs. 

The global heat-transfer measurements and predic- 
tions of heat transfer are presented in Table 4. 


Table 4. Measured and predicted heat transfer. 



Heat transfer, Btu/sec 

Experiments 

23.6 

Stratified flow analysis 

23.0 

Fully mixed analysis 

31.5 


The stratified-flow analysis predicts heat transfer rates 
that are within 2.6 % of experimental measurements. 
Given that the turbulence model was rather simplistic in 
nature and the turbulent Prandtl number was tailored to 
accurately reflect fiee-stream mixing conditions, a rather 
coarse, loosely packed grid was used, and the fact that the 
constituency of the stratified core itself was only modestly 
well defined, it is most likely that this high level of agree- 
ment between analysis and experiment is fortuitous. The 
predicted heat-transfer rates for the case of the fully-mixed 
awe are substantially higher than experimental measure- 
ments and that predicted in the stratified-flow analysis. 

Predicted and Measured Performa nce Parameters 

In the case of this particular film-cooled nozzle analy- 
sis, the prediction of nozzle performance was especially 
difficult since the chemical composition of the cote flow 
was not known and had to be assumed. Further, the mixing 
effects between the film and the fluid core were essential 
to predict specific impulse and characteristic velocity and 
would seem to depend heavily on the nature of the turbu- 
lence model. In light of this, the critical performance 
parameters of Ijp, c*, and Cp are compared for the analyti- 
cal and experimental results of the two nozzle cases. The 
results, shown in Table 5, indicate that the flow may have 
been even more stratified than what was assumed in the 


analysis as the stratified-flow predictions are somewhat 
higher than experimental measurements. The analysis pre- 
dicted a 6.4 % higher specific impulse and a 5.9 % higher 
characteristic velocity than was obtained experimentally. 
However, it should be noted that a different injector used 
in the test program yielded specific impulse and character- 
istic velocity measurements that were within 3 % of pre- 
dictions. Performance measurements for thrust coefficient 
are within 0.5 % of predictions for the stratified-flow anal- 
ysis case. 

In general, the stratified-core predictions for engine 
performance are about 5 % closer to the experimentally 
measured values than the fully-mixed core case. The 
exception is the nozzle vacuum thrust coefficient which is 
nearly the same in both analysis. Again, the accuracy of 
the Tethys model in predicting thrust coefficient and the 
relative insensitivity of this parameter with respect to flow 
stratification indicates that if the core flow had been more 
stratified than indicated in the analysis, it would have had 
little impact on the already very accurate thrust coefficient 
predictions. 

Table 5. Comparison of predicted and measured film- 
cooled nozzle performance parameters. 



Measured 

Fully mixed 

Stratified 

!sp 

304.1 

339.9 

323.5 

c* 

5947. 

6599. 

62%. 

C F 

1.645 

1.657 

1.653 


Transpiration Cooled Rocket Nozzle 

A hydrogen/oxygen transpiration-cooled rocket 
engine was tested at the NASA LeRC in the early 
1970’s. 20 The purpose of the program was to verify if 
hydrogen coolant could significantly reduce heat fluxes in 
the throat region of a rocket nozzle and to demonstrate if a 
new technique of creating a transpired region, via a series 
of fine platelet stacks, produced a reliable, efficient and 
economical means of transpiration cooling. Schematics of 
the baseline rocket engine and the transpiration-cooled 
spool piece section are shown in the Figures 18 and 19. 
The demonstration tests for hydrogen transpiration cooling 
were not performed on a conventional converging-diverg- 
ing rocket nozzle but rather with a plug-and-spool nozzle, 
with the transpiration coolant being injected into the main 
stream of the exhaust gases only along the spool piece and 
limited to a region 0.5 inches in length ending at the 
throat. The overall length of the spool piece was 4 inches 
while the gap between the spool piece and the plug was 
0.25 inches at the throat and 0.5 inches in the combustion 
chamber. The inside diameter of the spool piece was 2.6 
inches. 

The nominal operating pressure of the engine was 600 
psia and the nominal oxygen/fuel mixture ratio was 6. 
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Plug Spoolpiece 



Fig. 19. Transpiration cooled spoolpiece section. 

Several important assumptions were made with the 
Tethys code in modelling the flowfield of this engine. 
Because of the relatively incomplete nature of the experi- 
mental d at a reported, and the significant computational 
advantages, the hydrogen/oxygen propellant system was 
assumed to be comprised of a perfectly combusted mixture 
of steam and excess hydrogen. Furthermore, the tempera- 
ture of the transpired coolant was presumed to be equal to 
the wall temperature of the copper platelet stack through 
which it flowed. While no direct temperature measurement 
of the fluid being injected was made, a heat-transfer analy- 
sis performed during the experimental program indicated 
that the wall temperatures may have been 300 R hotter 
than the coolant fluid. However, an estimated temperature 
of the coolant fluid was used in the Tethys analysis to rep- 
resent the wall temperature in the transpired section, with 
the reasoning that the 300 R difference in temperature was 
far more critical to the change in enthalpy of the coolant 
than the conduction heat transfer. The temperature profile 
deduced from the experimental data and used in the Tethys 
analysis for the transpiration-cooled nozzle is shown in 
Figure 20. 



Fig. 20. Wall temperature distribution. 

An additional consideration investigated was the 
roughness effect of the transpired surface on the flow 
stream. Because the transpired region represented a series 
of rather coarse and rough (a surface roughness of 12 mils 
was measured) platelet stacks, a smooth- wall and a rough- 
wall analysis was performed to evaluate the effect of this 
condition. 

Mass Flux Predictions 

The mass flow rate required to maintain the wall at the 
prescribed temperature shown in Figure 20 is a critical 
parameter. Because the mass flow that enters a transpired 
region may exit via diffusion (either ordinary or thermal) 
and advection, the relative importance of each is also of 
interest In the figures that follow, the mass flux predic- 
tions for both the smooth-wall and the rough-wall analysis 
are presented. 

In Figure 21, the predictions of mass flux in the tran- 
spiration-cooled section of the nozzle are presented for the 
two analysis cases along with the averaged mass flow rate 
measured in the experimental program. One significant 
result is that the mass flow rates of a smooth, transpiration- 
cooled wall are appreciably lower than was experimentally 
observed. On the other hand, the rough-wall analysis 
appears to track the experimental mass fluxes rather well 
until reaching the center of the transpiration section. At the 
point where the prescribed wall temperature begins to 
increase (i.e., the center of the transpiration-cooled sec- 
tion), the mass injected into the fluid stream acts as a film 
coolant and serves to reduce the transpiration flow until 
reaching the edge of the transpiration-cooled region where 
the highly accelerating flow serves to dramatically 
increase coolant requirements, The extent to which this 
occurs may be overexaggerated because of the discontinu- 
ous nature of the adjoining boundary conditions. Similar 
trends are observed in the smooth-wall analysis but to a 
much lesser extent. Table 6 summarizes the differences 
between the observed and predicted transpiration mass 
flow rates in the plug-and-spool engine. Here it is 
observed that the smooth-wall prediction is nearly 1/2 that 
of experimental measurements, while the rough-wall pre- 
diction is much closer to experimental results, and is only 
13 % lower than experimentally measured flow rates. 
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Fig. 21. Predicted and experimental mass fluxes. 


Table 6. Total coolant flow rate requirements. 



Analysis 

Exp. 

Difference, % 

Smooth 

Rough 

— 

Flow rate, 
lbm/sec 

0.022 

0.039 

0.045 

Smooth: -51. 
Rough: -13. 


The individual contributions of mass flux into the fluid 
stream are shown in Figures 22 and 23 for the smooth- 
walled and rough-walled cases, respectively. 

Examination of the smooth-wall case shows many 
intgn-grin g features. The first is the importance of the Soret 
mass-flux term in both the solid-wall and transpiration- 
cooled sections of the nozzle. As the term is always posi- 
tive, the Soret term drives the lighter fluid (hydrogen) 
away from the wall and towards the mainstream of the 
mrhangt gase s, in the solid-wall regions of the flow, the 
positive Soret hydrogen flux term is balanced by the ordi- 
nary diffusion terms resulting in no net flux. On the other 
hand, in the transpiration-cooled section of the engine, the 
ordinary diffusion flux is in the positive direction with the 
net result that all hydrogen mass flux terms, advective and 
diffusive, act to drive the hydrogen into the flow stream. 



Fig. 22. Predicted hydrogen mass fluxes along the wall 
using a smooth-wall analysis. 


The rough-wall case exhibits many of the same trends 
as the smooth-wall case with the exception that in the tran- 


spiration-cooled section of the engine the contributions of 
hydrogen mass flux as the result of ordinary and thermal 
diffusio n are very small compared to the contributions of 
the advective and turbulent diffusion terms. This is to be 
expected as an extremely high level of surface roughness 
was used in the analysis. 



Fig. 23. Predicted hydrogen mass fluxes along the 
wall using a rough-wall analysis. 

The individual contributions of diffusion heat-flux 
components into the fluid stream are shown in Figure 24 
for the smooth-wall case. The results of the rough-wall 
case have qualitativ ely similar results. 5 The heat transfer 
rates due to conduction are seen to vary dramatically 
throughout the nozzle. The conduction heat-transfer rate 
decays significantly in the transpiration-cooled section of 
the nozzle, being about 20 Btu/(in. 2 -sec) just before the 
transpiration-cooled section and then decaying to less than 
10 Btu/(in. 2 -sec) in the transpiration-cooled section. The 
effect of this film cooling would probably be greater than 
suggested as the conduction heat flux would normally (i.e., 
if the wall were solid) increase by 50 % in this region as 
the result of flow acceleration 21 



Fig. 24. Predicted heat flux into the exhaust gases 
using a smooth-wall analysis. 

Unlike conduction, the flux of heat as the result of 
interdiffusional heat transfer is observed to be directed 
into the main stream of the flow. This is the result of the 
en thalp y of steam being negative at the prescribed wall 
temperatures. 

A final observation to be made is that the flux of heat 
as the result of Dufour heat transfer is negligible. 


10 




Conclusions 

An advanced engineering tool has been developed for 
the analysis of chemical hydrogen/oxygen rocket engines. 
Demonstration cases of a 1030: 1 area ratio nozzle, a film- 
cooled rocket engine, and a transpiration-cooled rocket 
engine, were analyzed. In the demonstration cases exam- 
ined, the model predicted thrust coefficient within 0.2 % to 

0. 5 . % of experimental measurements, which was well 
within experimental uncertainties of the test data. Heat- 
flux predictions also agree well with experimentally mea- 
sured values. Among some of its unique attributes, the 
model includes the effects of Soret mass transfer and 
Dufour heat transfer, where the Soret mass transfer term 
was shown to play an important role in the mass-flux con- 
tribution of a transpiration-cooled nozzle and in establish- 
ing the species mass fraction along the wall of a solid- 
walled engine. Conversely, the Dufour heat transfer term 
was shown to be negligible. 
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